PROCEEDINGS 

0F SCIENCE 




Critical examination of RHIC paradigms — mostly 
high p T 



o 

Michael J. Tannenbaum* 1 ^ 

Physics Dept., 510c, Brookhav en National Laboratory, Upton, NY 11973-5000,USA 
E-mail: mjt@bnl . gov 

CO 

A critical examination of RHIC paradigms is presented. Topics include: search for a critical point 
1) with a low energy scan; the lack of understanding of radiative processes in a medium in QCD 

compared in detail to examples from QED; the reason why some physicists started to measure 
^ particles at large pj in the 1960's; a review of the discovery of hard- scattering in p-p collisions in 

■ ^ i the 1970's via single-inclusive and two-particle correlations and application of these techniques at 

^1 RHIC. Several paradigms in both soft and hard physics which are popular at RHIC are discussed 

and challenged. 

CO 

m 

i— H 

00 

o 
o 



Workshop on Critical Examination of RHIC Paradigms 
April 14-17, 2010 
Austin, Texas, U.S.A. 



* Speaker. 

f Supported by the U.S. Department of Energy, Contract No. DE-AC02-98CH1-886. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos.sissa.it/ 



Critical examinations ofRHIC paradigms — mostly high pr 



Michael J. Tannenbaum 



1. The Quark Gluon Plasma (QGP) 

Let me start out with the famous schematic plots of the phase diagram of nuclear matter Fig. 1 [1, 
2]. The presumed trajectories of the evolution of the medium for collisions at RHIC and LHC cm. 



Figure 1: a) A schematic representation of the QCD Phase Diagram [1]. b) The location of the critical 
point, the separation between the lst-order transition and chemical freeze-out, and the focusing of the event 
trajectories towards the critical point, are not based on specific quantitative predictions, but are all chosen to 
illustrate plausible possibilities [2] 

energies are shown in Fig. la, where the axes are the temperature T vs. the baryon chemical poten- 
tial jig. The temperature for the transition from the Quark Gluon Plasma (QGP) to a hadron gas is 
taken as 170 MeV for jig = and the phase boundary is predicted to be a smooth crossover down to 
a critical point below which the phase boundary becomes a first order phase transition. Also shown 
are idealized trajectories for the RHIC cm. energy scan and future experiments at FAIR which are 
being performed in order to find the QCD critical point. 

One paradigm that I don't believe is that one can find the critical point by an energy scan at 
RHIC. Luciano Moretto at WWND2010 said, "critical fluctuations depend on a first order phase 
transition below the critical point whose large fluctuations drive the critical fluctuations", leading 
me to ask again whether all the interesting physics just outside the E802 aperture? What I mean 
by this is that since no fluctuations, other than from known correlations such as Bose-Einstein 
interference, were observed at the AGS fixed target Au+Au program at cm. energy ^/s^n = 4.8 
GeV (see below), or from the CERN Pb+Pb fixed target program at cm. energy y/sNN = 17.2 GeV 
(Fig. 2a) [3, 4] there is no first order phase transition to drive the critical fluctuations. Specifically, 
from AGS to RHIC energies, the non-random fluctuations of M PT , the event-by-event average pj, 
are about 1% of the (pt), averaged over all events, independent of y/sm (Fig- 2b) [5, 6] and show 
no evidence of fluctuations from a first-order phase transition. The only escape is if the QGP can 
not be reached in A+A collisions below the critical point, i.e. the endpoints of the trajectories in 
Fig. lb for y/sM = 17.2 GeV and below lie between the Chemical Freeze-out line and the phase- 
boundary. In principle, this is possible, but in my opinion it is very unlikely given the constancy of 
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Figure 2: a) Distribution of event-by-event average pr, M PT = pf=± YH=\ Pt { , where n is the number of 
particles on a given event, in 158 GeV/nucleon Pb+Pb collisions (data points); histogram is from mixed 
events [3]. Line is calculated random convolution of inclusive pr spectrum [4]. b) rms/mean of M PT distri- 
bution with random effect subtracted, Z DT = a / — ( ) as a function of x /smm [6] . 

V <Pr> \ {PT) Random 



Z pr as a function of ^/snn i n Fig- 2b. 



2. Given the Quark Gluon Plasma (QGP), what is my interest? 

The QGP is the only place in the universe where we can in principle and in practice understand 
QCD for color-charged quarks in a color-charged medium. How long will it take before we under- 
stand a quark in a QCD medium as well as we understand the passage of a muon through Copper 
in QED (Fig. 3)? 
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Figure 3: dE/dx of a /i + in Copper as a function of muon momentum [7]. 
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2.1 QED Bremsstrahlung in matter and crystals 

Even for such a seemingly simple reaction as muon bremsstrahlung, there is no simple for- 
mula [8]. For bremsstrahlung of an electron of energy E in a solid, apart from the standard Bethe- 
Heitler radiation (kdN/dk = d/X Q , where £=energy of radiated photon and d/X Q =the thickness in 
radiation lengths, X ), there are also very interesting effects [9]: the Landau-Pomeranchuk-Migdal 
(LPM) effect, which decreases the radiation at low x = k/E\ transition radiation at low x\ and a 
medium effect on the forward outgoing photon from Compton scattering off the atomic electrons 
which if inside the coherence length kills the dN/dk divergence. The coherence length is set by 
Wi- For bremsstrahlung of a photon with energy k by a particle with mass m, energy E: 



Lcoh 



— Rh — v ^min — 



m k 



2E(E-k) 



(2.1) 



where q = y |f | is the 4-momentum transfer to the target. Figure 4a [9] shows the bremsstrahlung 
spectrum of equivalent photons, kdN/dk, from 8 GeV electrons incident on a 3% radiation length 
(X ) Aluminum target. The pure Bethe-Heitler spectrum, kdN/dk = d/X Q is shown by the dashed 
line. Inclusion of the LPM suppression gives the dotted line. Further inclusion of the dielectric 



0.12 



a) 



0.1 



0.08 



o 
>< 



8> 0.06 - 



0.04 - 



0.02 



(c) 3% X Aluminum 

8 GeV beam 




0.2 



10 

k (MeV) 



1 "'■ ( 1 1 


"I 



















I 

5 

i 
I 

5 



"3 



100 500 



b) 



Figure 4: a) SLAC measurement of equivalent photon spectrum kdN/dk in 3% X Al target [9] illustrating 
LPM and dielectric suppression at low x = k/E. b) Calculation of coherent bremsstrahlung in a diamond 
crystal with a peak at 9.48 TeV for a proposed 10 TeV electron beam at the SSC [1 1]. 



suppression gives the histogram, which agrees with the measurement. For crystals, there is also 
coherent bremsstrahlung when the momentum transfer vector q equals a characteristic momentum 
of the reciprocal lattice (Bragg condition) which greatly increases the radiation [10] as shown in 
Fig. 4b [11]. 

I hope that the discussion in this section has illustrated the point that the understanding of 
dE/dx for quarks and gluons in a QGP, which is much more complicated than for leptons and 
photons in solids in QED, is in its infancy and has barely scratched the surface. 
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3. Why were some people studying "high px" physics in the 1960's 

The quick answer is that they were looking for a 'left handed' intermediate boson W ± , the 
proposed carrier of the weak interaction [12]. 

The first opportunity to study weak interactions at high energy was provided by the devel- 
opment of neutrino beams at the new accelerators in the early 1960's, the CERN-SPS and the 
BNL-AGS [13] where the /i-neutrino was discovered [14]. However, it was soon recognized that 
the intermediate (weak) boson, W ± might be more favorably produced in p-p collisions [15]. The 
canonical method for discovering the W-boson in p-p collisions was described by Nino Zichichi in 
a comment at the 1964 ICHEP, which deserves a verbatim quote because it was exactly how the 
W was discovered at CERN 19 years later [18, 19]: "We would observe the /i's from W-decays. 
By measuring the angular and momentum distribution at large angles of K and 7r's, we can predict 
the corresponding /i-spectrum. We then see if the /i's found at large angles agree with or exceed 
the expected numbers." I became aware of this "Zichichi signature" in my graduate student days 
when I checked the proceedings of the 12th ICHEP in Dubna, Russia in 1964 to see how my the- 
sis results were reported [16] and I found several interesting questions and comments by an "A. 
Zichichi" printed in the proceedings (Fig. 5). The W is still being measured by this method at 
RHIC (Fig. 6) [22] except now we use polarized p-p collisions and can measure the parity violation 



MUON-PROTON ELASTIC SCATTERING AT HIGH MOMENTUM TRANSFERS * 
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We have studied the elastic scattering 
of negative muons from liquid hydrogen at mo- 
mentum transfers of 550 MeV/c to 1050 MeV/c 
(q* = 7 to 26 fermi -3 }, using a detecting array 
of spark chambers and scintillation counters. 
The experiment was performed at the AGS acce- 
lerator of the Brookhaven National Laboratory, 
and the runs were divided into three stages, 
as shown below: 



of the proton in an aluminum plate spark cham- 
ber. One also measures the directions of the 
recoil proton and the recoil muon. This is ' 
equivalent to measuring three independent 
angles, from which one can infer for each event 
the value of k, and still overdetermine the scat- 
tering event by two degrees of freedom. This 
redundancy is used to select true scattering 
events from the background of false events, 
such as random coincidences, inelastic, u. — p 



It appears that, within the uncertainties 
of these preliminary results, the muon scatte- 
ring cross section for momentum transfers 
of up to I Gev/c is correctly described by the 
Rosenbluth formula and thee— p form factors. 
It is still too early to attempt a more quanti- 
tative definition of the possible deviation from 
the electron predictions. 



In connection with the problem of observing 
production of intermediate bosons, I would like 
mention that we have been studying at CERN t 
schemes: 



1) p+p-+W+ + W~ 



This process is described by the fallowing Feynman 

p/ W- 
Notice that this process is proportional to a 2 where 
a is the electromagnetic coupling constant. 

2) The second proposal studied would use the 
internal target of the proton synchrotron with 10 12 
protons per pulse incident onto the target. The pro- 
cess would be p + (J~^ -*-W* + anything. We wo- 
uld observe the fi'sfrom F-decays. By measuring the 



The accuracy of the highest point (statistical e: 
i i 1 1 i I r,o MeV/c, the e; 

a the point is about 15%. 



angular and momentum distribution at large angles 
of K and re's, we can predict the corresponding u.- 
spectrum. We then see if the u.'s found at large angles ■ 
agree with or exceed the expected number. A sup- 
plementary check can be made by measuring the 
polarization of these fi's. The polarization indicates 
the origin of these u.'s. Notice that the cross section 
for this process goes with where g is the f)-<lecay 
coupling constant. 



n t e 



uld like to use the fact that you are all tired 
in order to make a remark of linguistic rather than 
sicentific character. All the speakers used as nota- 
tions for neutral leptons the letters v„ and v^, This 
seems to be a very convenient notation. On the other 
hand, the terms which are usually used for neutral 
leptoiis-electron and muon neutrinos (and even «eie- 
ctron and muon type of neutrino!*), are too cumber- 
some. True, sometimes for v e the word ((neutrinos 
is used and for v^, the word «neutretto». The last 
term, however, is not very satisfactory since the last 
thirty years lost of particle ir.ehuiin:; strong inte- 
racting particles had been called that way. In addi- 
tion, it seems to me that both types of neutral lep- 
tons should conserve in their «name» the root «neu- 
trinos, which is widely associated with the unique 
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Figure 5: a)(left)Zichichi ICHEP 1964 [16, 17] and b)(right) W -> e + v from UA2 [20] 1982-1984 runs, 
about 1/4 of the final sample [21] 
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in W production which will enable us to measure flavor identified spin structure functions of u and 
d quarks [23]. Now back to the mid 1960's 
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Figure 6: PHENIX (a) and STAR (b) preliminary results [22] for W+ -> e+ + X from a short run at RHIC 
of polarized p-p collisions at = 500 GeV in 2009. The parity violation in real W + production has been 
measured for the first time with this data to ~ 4.0(7 significance, combined. 



3.1 The absence of high px single leptons leads to lepton-pair measurements. 

Proton-beam-dump experiments at the ANL-ZGS [24] and BNL-AGS [25, 26] looking for 
"large angle" muons didn't find any. Then the big question became, "How do you know how many 
W's should have been produced?". Chilton Saperstein and Shrauner [27] emphasized the need 
to know the timelike form factor of the proton in order to calculate the W production rate; and 
Yamaguchi [28] proposed that the timelike form factor could be found by measuring the number 
of lepton pairs (e + e~ or /i + /i~), "massive virtual photons", of the same invariant mass as the W 
(but also noted that the individual leptons from these electromagnetically produced pairs might 
mask the leptons from the W). This set off a spate of single and di-lepton experiments, notably the 
discovery by Lederman et al. of "Drell-Yan" pair production at the BNL-AGS [29, 30] followed by 
the experiments E70 and CCR at the new FNAL-Tevatron and CERN-ISR machines, respectively. 

The discovery of "Drell-Yan" pairs at the AGS proved to be seminal in future Relativistic 
Heavy Ion Physics as well as providing an interesting lesson. Figure 7a) shows the di-muon invari- 
ant mass spectrum da/dm^ from the collisions of 29.5 GeV protons in a thick Uranium target. 
There was definitely a dispute in the group about the meaning of the shoulder or "bump" or "??" 
for 2.5 < m m < 4.0 GeV/c 2 , which was apparently resolved adequately by the long forgotten the- 
ory paper [31] which produced the curve that beautifully agreed with the data (Fig. 7b) with the 
explanation/The origin of the shoulder comes from an interplay between the phase-space control 
of the integration region and the q 2 dependence of the coefficient." The important lesson that I 
learned from Fig. 7 is to NEVER be influenced by theoretical curves which "explain" your data. It 
is only when the curves fail to explain the data that you learn something definitive: the theory is 
wrong. 
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V I K W LETTERS 4 January 1971 




Figure 7: a) Dimuon invariant mass spectrum da/dm^^ [29]; b) theoretical prediction [31] for da/dm^^. 

Leon Lederman was very excited in 1970 to be in the possession of the di-muon continuum 
mass spectrum, da/dm^, because by combining this result at = 7.4 GeV with the newly found 
Bjorken Scaling [32] as used by Drell and Yan [30], he could calculate the W cross section at any 
y/s, and hence the sensitivities of his two proposals E70 [33] and CCR [34] (see Fig. 8). Details 
worthy of note from Fig. 8b are the e~ 6pT pion yield, the line with the worst imaginable background 
and the Jacobean peaks at pj = M w /2 from W — » e + X for various W masses. 



Proposal to The National Accelerator Laboratory 



E70-(F)NAL 



"Study of Lepton Pairs from Proton-Nuclear Interactions; 
Search for Intermediate Bosons and Lee-Wick Structure" 
W. LEE, L.M. LEDERMAN, J. APPEL, Columbia University, 
M. TANNENBAUM, Harvard University, L. READ, J. SCULLI, 
T. WHITE, and T. YAMANOUCHI, National Accelerator 
Laboratory. 

ABSTRACT 

We propose to observe lepton pairs emerging from hie 
energy proton-nuclear collisions. Large effective mass p 
probe the hadronic electromagnetic structure. The con- 
tinuum mass spectrum will be measured and any resonant 
structures in the mass range up to ~ 28 GeV will be 
detected with great sensitivity. The data provides a 
prediction, via Conserved Vector Current theory, for 
the production cross section for weak vector bosons and 
these are also sought in the mass range ~8-25 GeV. We 
also propose an initial photon-electron beam survey at 
high transverse momentum which is also a W-search with 
good sensitivity. 



June I?, 1970 + addendum Dec 1 970 -» 

Correspondent: L. M. Lederman, Columbia University 
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Figure 8: a) (left) Proposal for E70 Fermilab, June 17, 1970; b) (right) Cross section and background 
calculation for W ± — >> e ± + X from Addendum [33]. 
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4. The November Revolution and Birth of a QGP Paradigm 

Lederman's shoulder was explained in November 1974 by the discovery of the J/xy at the 
BNL-AGS [35] and at SLAC [36] (Fig. 9). This discovery revolutionized high energy physics since 
it was a heavy vector meson with a very narrow width (it decayed slowly) which implied a new 




Volume 33, Number 23 



PHYSICAL REVIEW LETTERS 



2 December 1974 



Experimental Observation of a Heavy Particle Jf 

J. J. Aubert, U. Becker, P. J. Biggs, J. Burger, M. Chen, G. Everhart, P. Goldhagen, 
J. Leong, T. McCorriston, T. G. Rhoades, M. Rohde, Samuel C. C. Ting, and Sau Lan Wu 

Laboratory for Nuclear Science and Department of Physics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 02139 



Y. Y. Lee 

Brookhaven National Laboratory, Upton, New York 11973 
(Received 12 November 1974) 

We report the observation of a heavy particle J, with mass m =3.1 GeV and width ap- 
proximately zero . The observation was made from the reaction £ + Be — e + + e~ + x by 
measuring the e *e~ mass spectrum with a precise pair spectrometer at the Brookhaven 
National Laboratory's 30-GeV alternating-gradient synchrotron. 



This experiment is part of a large program to 
study the behavior of timelike photons in p+p—e* 
+ e'+x reactions 1 and to search for new particles 
which decay into e*e~ and /i + m" pairs. 



daily with a thin Al foil. The beam spot size is 
3x6 mm 2 , and is monitored with closed-circuit 
television. Figure 1(a) shows the simplified side 
view of one arm of the spectrometer. The two 



Discovery of a Narrow Resonance in e + e Annihilation* 

J.-E. Augustin,t A. M. Boyarski, M. Breidenbach, F. Bulos, J. T. Dakin, G. J. Feldman, 
G. E. Fischer, D. Fryberger, G. Hanson, B. Jean-Marie, t R. R. Larsen, V. Liith, 
H. L. Lynch, D. Lyon, C. C. Morehouse, J. M. Paterson, M. L. Perl, 
B. Richter, P. Rapidis, R. F. Schwitters, W. M. Tanenbaum, 
and F. VannucciJ 

Stanford Linear Accelerator Center, Stanford University, Stanford, California 94305 
and 

G. S. Abrams, D. Briggs, W. Chinowsky, C. E. Friedberg, G. Goldhaber, R. J. Hollebeek, 
J. A. Kadyk, B. Lulu, F. Pierre, § G. H. Trilling, J. S. Whitaker, 
J. Wiss, and J. E. Zipse 

irence Berkeley Laboratory and Department of Physics, University of California, Berkeley, California 94720 
(Received 13 November 1974) 

We have observed a very sharp peak in the cross sec 
possibly (iV at a center-of-mass energy of 3.105±0.( 
full width at half-maximum is 1.3 MeV. 



Figure 9: (left) Sam Ting with / discovery at BNL; (right) / [35] and if/ [36] discoveries in PRL 

conservation law, similar to the discovery of strange particles in cosmic rays [37]. The J/\j/ was 
quickly understood to be a bound state of heavy c — c quarks (charmonium) [38] — the hydrogen 
atom of QCD. This was clear evidence for a "2nd" generation of quarks, and made all physicists 
believe in quarks and QCD. This discovery also changed the paradigm of di-lepton production 
from the measurement of background for W production to the search for new resonances. Greatly 
improved di-muon measurements by E70 in 1977 in p+A collisions at Fermilab [39] set the standard 
for di-muon production (Fig. 10a) and were rewarded with the discovery of an even heavier "3rd" 
generation of quark, the ^-quark [39]. At CERN, di-muon measurements were extended to h+A 
collisions by NA3 [40] and NA10 [41] and then to Pb+Pb collisions by NA50 (Fig. 10b) [42, 43]. 

4.1 J /y/ suppression as the first paradigm for detecting the QGP 

In 1986, Matsui and Satz [44] proposed that due to the Debye screening of the color potential 
in a QGP, charmonium production would be suppressed since the c — c quarks couldn't bind. This 
paradigm drove the field of RHI physics for two decades; and the observation of // \j/ suppression 
in A+A collisions [45, 43] is the CERN fixed target Heavy Ion program's claim to fame. How- 
ever, the interpretation is not straightforward because the // y is suppressed in p+A collisions (see 
Fig. 11). The PHENIX experiment at RHIC was specifically designed to measure J/ iff production 
at rest at mid-rapidity in p-p and Au+Au collisions [46] and found the same J /iff suppression at 
= 200 GeV as found at ^/s^ =17.2 GeV (158A GeV bombarding energy) at CERN by 
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a) ! ° ' b) M uu (GeV/c 2 ) 

Figure 10: a) CFS di-muon spectrum at ^/snn=27A GeV [39], with // y/, \j/ f , T clearly seen along with the 
Drell-Yan continuum, do /dm oc e~ lm0m . b) NA50 di-muon spectrum in 158A GeV Pb+Pb collisions [42] 

NA50 [43], further complicating the // ^-suppression 'paradigm'. Incredibly, the CERN fixed tar- 
get RHI program is still making progress, having (after 21 years [45]) finally measured the p+A 
cold nuclear matter effect for J/\j/ at the same incident energy, 158 GeV, as for Pb+Pb [47] (com- 
pare Fig. 11), which materially changes the result. At best, the the // ^-suppression paradigm is 
presently inconclusive. Hopefully, results from the LHC may settle the issue [48]. 




Figure 11: The J/y/ -+ jl+jl cross section in A+B collisions divided by A x B [42]. Line is (AB) a with 
a = 0.92 ± 0.01. Only the Pb+Pb data fall below this line. However, note the *. 
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5. A second paradigm for detecting the QGP, BDMPSZ [49, 50] 1997-1998 

In 1998, at the QCD workshop in Paris, Rolf Baier asked me whether jets could be measured in 
Au+Au collisions because he had a prediction of a QCD medium-effect (energy loss via gluon radi- 
ation induced by multiple scattering [51]) on color-charged partons traversing a hot-dense-medium 
composed of screened color-charges [52]. I told him [53] that there was a general consensus [54] 
that for Au+Au central collisions at ^Js^n = 200 GeV, leading particles are the only way to find jets, 



because in one unit of the nominal jet-finding cone, R = v /(At]) 2 + (A</>) 2 , there is an estimated 

The good news was that hard- scattering in p-p collisions had been discovered at the CERN 
ISR [55-57] by the method of leading particles, before the advent of QCD, and it was proved by 
single inclusive and two-particle correlation measurements in the period 1972-1978 that high pj 
particles are produced from states with two roughly back-to-back jets which are the result of scatter- 
ing of constituents of the nucleons as described by QCD, which was developed during this period. 
The other good news was that these techniques could be used to study hard- scattering and jets in 
Au+Au collisions and that the PHENIX detector had been designed to make such measurements 
and could trigger, measure and separate / and 7T° out to pj > 25 GeV in p-p and Au+Au collisions. 
In fact, in many talks, dating from as long ago as 1979 [58], I have been on record describing "How 
everything you want to know about jets can be found using 2-particle correlations". 1 

One lesson from this new paradigm is that a good probe of QCD in a fundamental system such 
as p-p collisions (Fig. 12) provides a well calibrated probe of QCD in more complicated collisions 
such as Au+Au. In the 1980's when RHIC was proposed, hard processes were not expected to play 
a major role in A+A collisions. However, starting in 1997 [61], inspired by Rolf and collaborators, 
and before them by the work of Gyulassy [62] and Wang [63], I indicated that my best bet on 
discovering the QGP was to utilize semi-inclusive 7T° or production in search for "high pj 
suppression". 

5.1 Soft vs. hard-physics in p-p collisions 

Fig. 12a nicely illustrates two important points for 7T° production in p-p collisions. Apart from 
the beautiful data, the main emphasis is the agreement of NLO pQCD with the measurement for 
/?r>3 GeV/c. Actually this surprised some of my colleagues from the lepton- scattering commu- 
nity, but was no surprise to me because such data in the 1970's is what made people believe in QCD 
(see Sec. 7 below). The other important issue is the transition from the exponential "soft physics" 
spectrum for pj < 1 .5 GeV/c to the power law spectrum of hard scattering at larger pj. The e~ 5 - 5pT 
exponential at y/s = 200 GeV, is very close to the original cosmic-ray "Cocconi Formula" [66], 
e~ 6pT , used by Lederman (Fig. 8b) in 1970, but the power-law hard-scattering spectrum for p T > 3 
GeV/c is much larger than the "worst imaginable background". This is the hard- scattering com- 
ponent which I shall discuss again below but I first wish to emphasize some important issues in 
Relativistic Heavy Ion physics which are dominated by the soft physics. 

Recently, I had to amend this statement to 'almost everything' because we found that the fragmentation function 
can not be measured in di-hadron correlations because the two particles, e.g. 7T° — h, are both fragments of jets [59, 60] 
(see Sees. 7.3, 8.2.1 below). 
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Figure 12: (left) PHENIX measurement of invariant cross section, Ed 3 a/d 3 p, as a function of transverse 
momentum pr for 7T° production at mid-rapidity in p-p collisions at cm. energy y/s = 200 GeV [64]. (right) 
PHENIX measurement of 7T° in p-p collisions at y/s = 62.4 GeV [65]. 

6. The important of soft-physics in Relativistic Heavy Ion Collisions 

In dealing with Relativistic Heavy Ions, there are three dramatic differences from p-p physics. 




Figure 13: a) (left) A p-p collision in the STAR detector viewed along the collision axis; b) (center) Au+Au 
central collision at y/sNN = 200 GeV in the STAR detector; c) (right) Au+Au central collision at ^/snn = 200 
GeV in the PHENIX detector. 

• The particle multiplicity is ^ A times larger in A+A central collisions than in p-p collisions 
(Fig. 13). This is why there is such a huge energy, ~ 375 R 2 GeV, in a jet cone of radius R at 
y/s^ = 200 GeV. 
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• There are huge azimuthal anisotropies in A+A which don't exist in p-p which are interesting 
in themselves, and are useful, but sometimes can be troublesome. 

• Space-time issues, both in momentum space and coordinate space are important in RHI. For 
example: 

1. When in time and space does a parton fragment? Is this different for light and heavy 
quarks? When are particles formed? 

2. The Dokshitzer textbook formula, Tp = ER 2 , should, I think, be corrected to Tp = 
yR = RxE /M = ERXc, where R, M, are the radius and mass of the particle and Ac its 
Compton wavelength. 

3. Would a proton embedded in a QGP dissolve? How long does this take? How is this 
related toJ/if/ suppression? 

6.1 Soft physics dominates multi-particle production in both p-p and A+A collisions 

A quantity related to multiplicity, but more relevant for the study of jets, is the distribution of 
transverse energy, E T = ^Z^sinfy, where the sum is taken over all particles emitted on an event 
into a fixed but large solid angle, typically measured in a calorimeter, corrected for the calorimeter 
response and then corrected to the solid angle At] = 1, A0 = 2k [67]. Ej is also related to M PT , 
the event-by-event average pj (Fig. 2a), except that it is the sum of pj of all particles on every 
event rather than the average. dE T /dy is thought to be related to the co-moving energy density in 
a longitudinal expansion [69, 68], and taken by all experiments as a measure of the energy density 
in space, "the Bjorken energy density", 

£Bj - dy T F 7TRi (6A) 

where T/7, the formation time, is usually taken as 1 fm/c, 7iR 2 = A± is the transverse overlap area 
of the collision, and d (E T ) /dy is the co-moving energy density. 

The main importance of Ej distributions in RHI collisions (see Fig. 14a) is that they are 
sensitive primarily to the nuclear geometry of the reaction, and hence can be used to measure 
the centrality of individual interactions on an event-by-event basis. 2 Thus, the transverse energy 
or multiplicity distribution in A+A collisions has a characteristic shape which is quite different 
from that in p-p collisions since it is sensitive to the number of participants or centrality of an 
A+A collision and integrates over all impact parameters. This is shown in Fig. 14a for Au+Au 
collisions at AGS (y/sm = 5.4 GeV) and RHIC (^/s^ = 200 GeV), together with an estimate 
at RHIC for p-p collisions. The fact that the scaled AGS and RHIC Ej distributions lie one on 
top of each other shows that the shape of E T distributions is essentially entirely dominated by the 
nuclear-geometry at both the AGS and RHIC; which implies that there are no critical fluctuations 
at either the AGS or RHIC. Above the upper 0.5 percentile of the distribution (E T ^710 GeV on 
Fig. 14), measurements in smaller solid angles show increasingly larger fluctuations represented by 
flatter slopes above the knee of the distribution; but this is a random fluctuation (fewer particles in 



2 However, in PHENIX, centrality is determined far from mid-rapidity to avoid possibly biasing any of the mid- 
rapidity measurements. 



12 



Critical examinations ofRHIC paradigms — mostly high pr 



Michael J. Tannenbaum 



PHENIX and E802 E T Transverse Energy corr to A7j = l and A*=27T 



a) 




Figure 14: a) E T distributions in Au+Au collisions: PHENIX at y/sm = 200 GeV [70] compared to the 
E802 (non-response-corrected) Ej at ^/snn = 5.4 GeV [71, 72] scaled up by 8.1. The quantity plotted 
is the normalized yield, dY/dEr (GeV -1 ), which integrates to 1, i.e J dEj dY jdEj = 1. Different filled 
symbols represent PHENIX measurements in solid angles At] = 0.76, A(f) = m x tt/8 at mid rapidity with 
m— 1, ...5. Dashed line is an estimate of Ej distribution in p-p collisions at yj~s = 200 GeV for m = 5. b) 
PHENIX m = 5 raw Ej distribution, EjEMC, at s/snn = 200 GeV together with calculated distributions for 
participant scaling, WNM (Et Af part /2), and Ncoll scaling (Et N co \i) based on convolutions of the p-p 
distribution. 

the smaller solid angle fluctuate more). This is visible only for the most central collisions when the 
nuclear-geometrical fluctuations have saturated (all nucleons participating). The Bjorken energy 
density can be straightforwardly computed from from Fig. 14a at RHIC (^s^n = 200 GeV), with 
£Bj x t f = 5.4 ±0.6 GeV/fm 2 [70] for the 5% most central collisions (E T > 580 GeV). For the 
AGS data, the ratio of 8.1 shown in Fig. 14 can not be used to get £#/ x x F at y^Aw = 5.4 GeV 
because the E802 data have not been corrected for the calorimeter response, which is much more 
difficult in the fixed target geometry. However, Ref. [70] calculated the fully corrected E802 E T 
from charged particle measurements and obtained Egj x ?f = 1.0 GeV/fm 2 at ^/snn = 4.8 GeV. 

6.2 Another paradigm we should do away with 

It has been popular to characterize the average values of dE T /dr] (or the charged multiplicity 
density dn ch /dr\) in A+A collisions as a function of centrality by the equation [73, 74] 

dnfa/dr} = (1 -x) dn%/dr\ x Af par t/2 + x dnf p /dr\ x 7V coll . (6.2) 

While this may seem reasonable for the average value (drf^/di]) at a given centrality, is nonsense 
for the distribution, which if such a representation were true would be the weighted sum of (1 — x) x 
the WNM curve + x x the Ncoll curve in Fig. 14b, which obviously looks nothing like the actual 
Ej EMC measurement. (A more reasonable representation is dn^/dr] oc N part a .) Also the claim 
that the Ncoll scaling of the soft p-p E T or multiplicity n C p p distributions (or average values) in 
Eq. 6.2 has anything to do with hard- scattering [74] is equally nonsensical. This could be further 
elaborated by looking for jets in the p-p data; but such measurements (not finding jet activity until 
4-5 orders of magnitude down in Ej cross section) were already done by COR [75] and UA2 [76] 
25 years ago (see Fig. 21 below). 
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6.3 Another soft-physics paradigm to test at the LHC — Hydrodynamics 

In my opinion, the most interesting soft physics issue for the LHC concerns the possible in- 
crease of the anisotropic flow V2 beyond the 'hydrodynamic limit'. Wit Busza's extrapolation [77] 
of V2 to the LHC energy is shown in Fig. 15a, a factor of 1.6 increase from RHIC. A previous 
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Figure 15: a) (left) Busza's extrapolation of V2 to LHC [77]. b) (center) v^jz vs 'Bjorken multiplicity 
density', (l/S)dN c \ l /dy [78]. c) (right) 'Hydro Limit' calculated in viscous Hydrodynamics for several 
values of the initial energy density eo [79]. 

paper by NA49 [78] which compared V2 measurements from AGS and CERN fixed target exper- 
iments to RHIC as a function of the 'Bjorken multiplicity density', dn^/di] /S, where S = is the 
overlap area of the collision zone, showed an increase in V2/£ from fixed target energies to RHIC 
leading to a "hydro limit", where £ is the eccentricity of the collision zone (Fig. 15b). This limit 
was confirmed in a recent calculation by Uli Heinz using viscous relativistic hydrodynamics [79] 
which showed a clear hydro-limit of V2/£ = 0.20 (Fig. 15c). This limit is sensitive to the ratio of 
the viscosity/entropy density, the now famous r]/s, but negligibly sensitive to the maximum energy 
density of the collision, so I assume that this calculation would give a hydro-limit at the LHC not 
too different from RHIC, V2/£ ~ 0.20. Busza's extrapolation of a factor of 1.6 increase in V2 from 
RHIC to LHC combined with V2/£ from Fig. 15b gives V2/£ = 0.32 at LHC. In my opinion this 
is a measurement which can be done to high precision on the first day of Pb+Pb collisions at the 
LHC, since it is high rate and needs no p-p comparison data. Personally, I wonder what the hydro 
aficionados would say if both Heinz' and Busza's predictions were correct? 

7. Hard-Scattering in p-p collisions 

We now go back to where we were in the search for high p T leptons and di-leptons before we 
got diverted by the // \jf discovery (and QGP paradigm) and soft-physics issues. 

7.1 Bjorken scaling and the parton model 

The idea of hard- scattering in p-p collisions dates from the first indication of pointlike structure 
inside the proton, in 1968, found at SLAC by deeply inelastic electron-proton scattering [80], i.e. 
scattering with large values of 4-momentum transfer squared, Q 2 , and energy loss, V. The discovery 
that the Deeply Inelastic Scattering (DIS) structure function 

F 2 (Q 2 ,v)=F 2 (^ r ) (7.1) 
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"scaled" i.e just depended on the ratio 




independently of Q 2 , as originally suggested by Bjorken [32], led to the concept of a proton com- 
posed of point-like 'partons' [81]. The deeply inelastic scattering of an electron from a proton is 
simply quasi-elastic scattering of the electron from point-like partons of effective mass Mx, with 
quasi-elastic energy loss, V = Q 2 /2Mx (Fig. 16a). The probability for a parton to carry a fraction x 
of the proton's momentum is measured by ^(x) jx. 




Figure 16: a) Schematic e-p DIS via e-parton scattering [81] b) Schematic parton-parton scattering in p-p 
c) Predicted cross section according to reaction in (b). 

Since the partons of DIS are charged, and hence must scatter electromagnetically from each 
other in a p-p collision (Fig. 16b), Berman, Bjorken and Kogut (BBK) [82] and subsequent au- 
thors [83, 84], derived a general formula for the cross section of the inclusive reaction 

p + p^C + X (7.3) 

using the principle of factorization of the reaction into parton distribution functions for the protons, 
fragmentation functions to particle C for the jet of particles from the scattered parton and a short- 
distance parton-parton hard scattering cross section. The invariant cross section for the inclusive 
reaction (Eq. 7.3), where particle C has transverse momentum p T near mid-rapidity was given by 
the general 'scaling' form: [83] 

£ ^ = ^ F (^) = ^F G (^) where xt = 2 Pt /^s . (7.4) 

The cross section has 2 factors, a function F (G) which 'scales', i.e. depends only on the ratio of 
momenta; and a dimensioned factor, Pj Uefi (y/~s~ n&fi ), where /z e ff gives the form of the force-law 
between constituents. For QED or Vector Gluon exchange [82], n e ff = 4, and for the case of quark- 
meson scattering by the exchange of a quark [83], ft e ff-8. Using this formalism, BBK predicted 
(Fig. 16c) hard- scattering must exist in p-p collisions since the charged partons of DIS must scatter 
electromagnetically, this "electromagnetic contributions may be viewed as a lower bound on the 
real cross section at large pj\ 
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7.2 Discovery of high pr n production in p-p collisions 

The CERN Columbia Rockefeller (CCR) Collaboration [55] (and also the Saclay Strasbourg [56] 
and British Scandinavian [57] collaborations) at the CERN-ISR measured pion production over a 
large range of transverse momenta, unavailable in cosmic ray studies or at lower y/s. The e~ 6pT 
soft- spectrum at low pj which depends very little on y/s breaks to a power law at high p T with a 
strong and characteristic y/s dependence (Fig. 17a). The large rate at high p T , much larger than 
expected in Figs. 8b and 16c, indicates that partons interact strongly EM) with each other. 
However, the QED or Vector Gluon form of BBK [82], pj 4 F(pr/y/s), was not observed in the 
experiment but scaling with ft e ff — 8 was observed (Fig. 17b), which spawned a whole host of new 
theories. [83]. The first application of QCD to wide angle hadron collisions [84], showed that pure 
xt scaling breaks down: n e ff varies according to the xj and y/s regions used in the comparison, 
^eff —> n G ff(xT,y/s); but this was not enough to resolve n e ff — 8 issue, which turned out to be due to 
"pr broadening" by initial state parton transverse momentum, the "kj effect." [85]. 




Pl /Vs 



Figure 17: a) (left) CCR [55] transverse momentum dependence of the invariant cross section at five center 
of mass energies, b) (right) The same data multiplied by p\, using the best fit value of n = 8.24 ± 0.05, with 
F = Ae~ bx± , plotted vs p± / y/s. 

7.3 Everything you want to know about jets can be found using 2-particle correlations. 

Following the discovery of hard- scattering in p-p collisions at the CERN-ISR by single parti- 
cle inclusive measurements, the attention of experimenters turned to measuring the predicted di-jet 
structure of the hard-scattering events using two-particle correlations. At the CERN-ISR, from 
1975-1982, two-particle correlations showed unambiguously that high pj particles come from di- 
jets. The outgoing jet-pair of hard-scattering obeys the kinematics of elastic- scattering (of partons) 
in a parton-parton cm. frame which is longitudinally moving with rapidity y = (l/2)ln(jti/jt2) 
in the p-p cm. frame. Hence, the jet-pair formed from the scattered partons should be co-planar 
with the beam axis, with two jets of equal and opposite transverse momentum. Thus, the outgoing 
jet-pair should be back-to-back in azimuthal projection. It is not necessary to fully reconstruct the 
jets in order to measure their properties. In many cases two-particle correlations are sufficient to 
measure the desired properties. Many ISR experiments provided excellent 2-particle correlation 
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measurements [86]. However, the CCOR experiment [87] was the first to provide charged parti- 
cle measurement with full and uniform acceptance over the entire azimuth, with pseudorapidity 
coverage — 0.7 < ?] < 0.7, so that the jet structure of high p T scattering could be easily seen and 
measured. 



0.3<p T <1.0 0.3<p T <1.0 




A(|> A(f> - Xe 



Figure 18: a) (left) trigger-side (A<j> = ±n/2) b) (left-center) away-side (A(j> = ±7r/2 around 7t rad.) cor- 
relations of charged particles with indicated pr = pr a for 7t° triggers with prt > 7 GeV/c, for 5 intervals of 
Pt ; c) (right-center) xe distributions from this data, d) right) [top] Comparison [90] of away side charged 
hadron distribution triggered by a 7T° or a jet, where z K o = xe and Zj = PT a /pT a - [bottom] same distributions 
with 7T° plotted vs z!j = (zt)x£. 

The azimuthal distributions of associated charged particles relative to a 7T° trigger with trans- 
verse momentum p Tt > 1 GeV/c are shown (Fig. 18a,b) for five intervals of associated particle 
transverse momentum p T . In all cases, strong correlation peaks on flat backgrounds are clearly vis- 
ible, indicating the di-jet structure which is contained in an interval A0 = ±60° about a direction 
towards and opposite to the trigger for all values of associated pj (> 0.3 GeV/c) shown. The small 
variation of the widths of the away-side peaks for p T > 1 GeV/c (Fig. 18b) indicates out-of-plane 
activity of the di-jet system beyond simple jet fragmentation (the kj effect [85]). 

Following the methods of previous CERN-ISR experiments [88, 89] and the best theoretical 
guidance [90], the away jet azimuthal angular distributions of Fig. 18b, which were thought to 
be unbiased, were analyzed in terms of the two variables: p out = /?rsin(A0), the out-of-plane 
transverse momentum of a track, and 

-p T -pT t -Pt cos(A0) z irT j T ,r^!!^ L Wi^JET 

\PT t r PT t Zt ^J^_ _Jk T 

where it — PT t /pT t is the fragmentation variable of the trigger jet with pj t , and z a — PT a /pT a is 
the fragmentation variable of the away jet (with pj a ). Note that xe would equal the fragmenation 
fraction z of the away jet, for z t -> 1, if the trigger and away jets balanced transverse momentum, 
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i.e. if jth = PT a /pT t = 1. It was generally assumed, following the seminal article of Feynman, Field 
and Fox [90], that the pj a distribution of away side hadrons from a single particle trigger [with pj t \ 
corrected for (z t ), would be the same as that from a jet-trigger and follow the same fragmentation 
function as observed in e + e~ or DIS [91]. The xe distributions [87, 92] for the data of Fig. 18b 
are shown in Fig. 18c and show the fragmentation behavior expected at the time, e~ 6z ~ e - 6x E(zt) > 
Also, relevant to a recent claim [93], Fig. 18c from 1979 showed that there were no di-jets each of 
a single particle, as claimed by another ISR experiment [94], since there is no peak at xe = 1. 

As noted previously 1 , the title of this section should start with "Almost" because we subse- 
quently learned at RHIC [59] that the xe distribution from di-hadrons (Fig 18c) does not measure 
the fragmentation function (see Sec. 8.2.1 below). 

7.4 Direct searches for Jets 1975-1982: false claims, confusion, skepticism; finally success. 

One of the original reasons for measuring jets in p-p collisions was that it was assumed that 
the jet cross section would measure the parton cross section without reference to the fragmentation 
functions. However, due to the fact (which was unknown in the 1970's) that jets in 4tt calorimeters 
at ISR energies or lower are invisible below y/§ ~ E T < 25 GeV [100], there were many false claims 
of jet observation in the period 1977-1982 [101]. The coup-de-grace to the concept of finding jet 
structures in large aperture calorimeters was provided by the first measurement of an E T distribu- 
tion in the present usage of the terminology (see Sec. 6.1 above) by the NA5 experiment at CERN 
(Fig. 19a) [102] who found that "The events selected by the full calorimeter trigger show no dom- 
inant jet structure. They appear to originate from processes other than two constituent scattering." 
This led to mass confusion in the hard-scattering community in the period 1980-1982 and to great 




EJGeV) 

Figure 19: a)(left) NA5 structureless E T distribution [102]. b,c) UA2 jet event from 1982 ICHEP [103]: 
b)(center) event shown in geometry of detector; c)(right) "Lego" plot of energy in calorimeter cell as a 
function of angular position of cell in polar (6) and azimuthal (<E>) angle space. 

skepticism about jets in hadron collisions, particularly in the USA. A 'phase change' in belief- 
in-jets was produced by one UA2 event at the 1982 ICHEP in Paris [103], which, together with 
the first direct measurement of the QCD constituent- scattering angular distribution, E^(cos0*) 
(Eq. 7.5), using two-particle correlations, presented at the same meeting (Fig. 20) [103, 104], gave 
universal credibility to the pQCD description of high pj hadron physics [105, 91, 106]. Since that 
time, QCD and jets have become the standard tools of High Energy Physics. 
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Figure 20: a) (left 3 panels) CCOR measurement [103, 104] of polar angular distributions of 7T° pairs with 
net pr < 1 GeV/c at mid-rapidity in p-p collisions with = 62.4 GeV for 3 different values of % 71 invariant 
mass M KK . b) (rightmost panel) QCD predictions for Z^(cos 0*) for the elastic scattering of gg, qg, qq' , qq, 
and qq with cc s (Q 2 ) evolution. 

As a final note, recall that hard-scattering is visible by a break in the exponential single particle 
inclusive spectrum after roughly 3 orders of magnitude in cross section (Fig. 12). However, finding 
jet structure in Ej distributions measured in An calorimeters is much more difficult because the 
jets are only visible by a break in the E T spectrum after 5-7 orders of magnitude of cross section 
(Fig. 21) [76, 75], since the Ej spectrum is even more dominated by soft physics than the single 
particle spectrum. 




Fitjure 2 Distribution of the total transverse energy V J, T observed in l.he UA2 central 

calorimeter. E T (GeV) 

Figure 21: E T measurements by: a)(left) UA2 in p-p collisions at = 540 GeV [76]; b)(right) CCOR in 
p-p collisions at ^Js = 62.3 GeV [75]. The lines indicate extrapolation of the nearly exponential spectrum at 
lower Ej . 
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7.5 Status of theory and experiment, circa 1982 

Hard- scattering was visible both at ISR and FNAL (Fixed Target) energies via inclusive single 
particle production at large pj > 2-3 GeV/c. Scaling and dimensional arguments for plotting data 
revealed the systematics and underlying physics. The theorists had the basic underlying physics 
correct; but many important details remained to be worked out, several by experiment. The trans- 
verse momentum imbalance of outgoing parton-pairs, the "/^-effect", was discovered by experi- 
ment [89, 95], and clarified by Feynman, Field and Fox [90]. The first modern QCD calculation 
and prediction for high p T single particle inclusive cross sections, including non-scaling and the kj 
effect was done in 1978, by Jeff Owens and collaborators [96] closely followed and corroborated 
by Feynman and collaborators [97], under the assumption that high p T particles are produced from 
states with two roughly back-to-back jets which are the result of scattering of constituents of the 
nucleons (partons). The overall p + p hard- scattering cross section in "leading logarithm" pQCD 
is the sum over parton reactions a + b -» c + d (e.g. g + q — > g + q) at parton-parton center-of-mass 
(cm.) energy \f§ = yJx\X2S. 

, f° „, = =~ EAWAW ^£»»(cosr) (7.5) 

dx\dx2dcos0* dsdydcosQ* 2x\X2 

where f a (x\), fbfa), are parton distribution functions, the differential probabilities for partons a 
and b to carry momentum fractions x\ and X2 of their respective protons (e.g. ufa)), and where 0* 
is the scattering angle in the parton-parton cm. system. The parton-parton cm. energy squared is 
s = x\X2S, where is the cm. energy of the p-p collision. The parton-parton cm. system moves 
with rapidity y = (l/2)ln(jci/jC2) in the p-p cm. system. Equation 7.5 gives the pj spectrum of 
outgoing parton c, which then fragments into a jet of hadrons, including e.g. 7T°. The fragmentation 
function D^°(z) is the probability for a 7T° to carry a fraction z = p K ° /p c of the momentum of 
outgoing parton c. Equation 7.5 must be summed over all subprocesses leading to a 7C° in the 
final state weighted by their respective fragmentation functions. In this formulation, f a (xi), fb{x2) 
and tif \z) represent the "long-distance phenomena" to be determined by experiment; while the 
characteristic subprocess angular distributions, E^(cos0*), and the coupling constant, oc s (Q 2 ) = 
^ln(Q 2 /A 2 ), are fundamental predictions of QCD [98, 99]. 



8. Hard-Scattering at RHIC 

For the past decade the single and two-particle techniques developed in the period 1972-1982, 
largely at the CERN-ISR were used exclusively at RHIC to study hard-scattering in p-p A+A 
collisions, with outstanding results. 

8.1 Single particle (semi-) inclusive measurements 

The p-p single particle inclusive pr spectra measurements in Fig. 22 clearly illustrate the im- 
portance of particle identification since all the different particles have different systematic behavior. 
PHENIX concentrates on mesons and photons since they provide the main backgrounds to lepton 
and photon searches for which we built this specialized detector, while STAR concentrates on hy- 
perons. Also, inclusive single particle measurements have very precise ~ 1% pj scales in contrast 
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to jets. As shown in Fig. 22b for p-p collisions, jets are excellent if you can do it because they give 
the best rates. My favorite jet measurement is a search for parity violation in 100 GeV jets to look 
for quark substructure [107]. Fig. 22a, on a semi-log scale, emphasizes measurements at the lowest 
Pt including the soft-physics range and the transition to hard-scattering, while the log-log scale on 
Fig. 22b emphasizes the power-law behavior at high p T . It is interesting to note that the 71° and 
jet spectra are parallel to each other (the same p\ dependence, see Eq. 8.3 below), while the direct 
7 are flatter, corresponding to an n e ff closer to 4 on an xj scaling plot. There should be a simple 
explanation, but I don't know it. 




Figure 22: PHENIX pr invariant cross section at mid-rapidity in p-p collisions at ^/s = 200 GeV. a) Com- 
pilation of mesons needed for understanding leptonic background [108] b) Direct-/ [109, 110],7T° [64] and 
STAR Jets [111] for hard- scattering studies. 



8.1.1 High p T suppression in Au+Au collisions 

On of the major, if not the major discovery at RHIC is the huge suppression of particle produc- 
tion at large pj in A+A collisions. So far this appears to be consistent with the QCD predictions 
of energy loss or absorption of the outgoing hard- scattered color-charged partons due to interac- 
tions with the presumably deconfined and thus color-charged medium (a liquid version of the QGP) 
produced in central A+A collisions at RHIC [112]. The suppression is represented by the nuclear 
modification factor, the ratio Raa of the measured semi-inclusive yield (e.g of 71) for a given cen- 
trality class to the p-p cross section scaled by (Taa) the average overlap integral of the nuclear 
thickness functions for that centrality class (Eq. 8.1). For pure point-like hard-scattering, Raa = 1- 

= d^/pTdprdyNZ 1 
M {T AB ) x d 2 o* /p T dp T dy 
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The striking differences of Raa(pt) in central Au+Au collisions for the many particles measured 
by PHENIX at RHIC (Fig. 23) illustrates the importance of particle identification for understanding 
the physics of the medium produced at RHIC. For the entire region p T < 20 GeV/c so far measured 




16 18 

P T (GeV/c) 



Figure 23: Raa compilation from PHENIX for central Au+Au collisions. With the exception of the direct- 7 
internal-conversion data (7*) for pj < 4 GeV/c where a fit to the p-p data (which is not exponential) is used, 
all the other values of Raa are computed from the measured Au+Au and p-p pr spectra. 

at RHIC, all particles show suppression for pj > 3 GeV/c with the exception of p + p, which are 
enhanced in the region 2 < p T < 4 GeV/c (the "baryon anomaly" [113]), and the direct 7, which 
are not suppressed, presumably because the outgoing 7's do not interact with the color charged 
medium. New and truly notably this year is the behavior of Raa of direct- 7 for p T < 2 GeV/c [109], 
which is totally and dramatically different from all other particles, exhibiting an order of magnitude 
exponential enhancement as pj 0. This exponential enhancement is certainly suggestive of a 
new production mechanism in central Au+Au collisions different from the conventional soft and 
hard particle production processes in p-p collisions and its unique behavior is attributed to thermal 
photon production by many authors (e.g. see citations in reference [1 14]). 

8.2 Two-particle correlations at RHIC 

8.2.1 A RHIC discovery-di-hadron correlations do not measure the fragmentation function 

The two-particle (di-hadron) correlations as measured at RHIC follow the systematics from 
previous measurements at the CERN-ISR as already shown in Sec. 7.3. However, we learned some- 
thing new and important at RHIC [59]-the di-hadron correlation, e.g. 7T -h, where both hadrons 
are fragments of jets, does not measure the fragmentation function. 

PHENIX [59] attempted to measure the mean net transverse momentum of the di-jet, (pr pa i r ) = 
y/2(kr), in p-p collisions at RHIC, where kj represents the out-of-plane activity of the hard- 
scattering [90, 91]. This requires the knowledge of (z t ) of the trigger 7T°, which PHENIX attempted 
to calculate using a fragmentation function derived from the measured xe distributions (Fig. 24a). 
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Figure 24: a) xe distributions from PHENIX [59] for 7T°-h correlations in p-p collisions at ^fs = 200 GeV 
for several values of pj t . The solid and dashed lines represent calculations of the distribution for quark 
(solid lines) and gluon (dashed lines) fragmentation functions based on exponential fits to the LEP measure- 
ments [115, 116] shown in b). See Ref. [59] for details 



It didn't work. Finally, it was found that starting with either the quark « exp(— 8.2 • z) or 
the gluon « exp(— 11.4 • z) fragmentation functions from LEP (Fig. 24b solid and dotted lines), 
which are quite different in shape, the results obtained for the xe distributions (solid and dotted 
lines on Fig. 24a) do not differ significantly! Although nobody had noticed this for nearly 30 
years, the reason turned out to be quite simple. The integration over z t of the trigger jet for fixed 
trigger particle pr t is actually an integral over the trigger jet pj t . However since the trigger and 
away-jets are always roughly equal and opposite in transverse momentum, integrating over pj t 
simultaneously integrates over p Ta and thus also integrates over z of the away -jet. 

With no assumptions other than a power law for the jet pj t distribution (do q / pT t dpT t =Ap^ n ) , 
an exponential fragmentation function (D^(z) = Be~ bz ), and constant x^, for fixed pj t as a function 
of pr a , it was possible to derive the xe distribution in the collinear limit, where pr a = %EPT t [59]: 



dPn 

dxE 



&N(n-l)^r - * , (8.2) 



and N = (m) is the multiplicity of the unbiased away -jet. The shape of the xe distribution (Fig. 24a) 
is given by the power n of the partonic and inclusive single particle invariant transverse momentum 
spectra and does not depend on the exponential slope of the fragmentation function (Fig. 24b). Note 
that Eq. 8.2 provides a relationship between the ratio of the away and trigger particle's transverse 
momenta, xe ~ PT a /pr t , which is measured, to the ratio of the transverse momenta of the away 
and trigger jets, xh = pr a / Pr t , which can thus be deduced. In p-p collisions the imbalance of the 
away-jet and the trigger jet (xh ~ 0.7 — 0.8) is caused by /^-smearing [90, 59]. In A+A collisions, 
Xh is sensitive to the relative energy loss of the trigger and associated jets in the medium, which can 
thus be measured [85]. 

The same derivation gives a simple formula for Bjorken's parent-child ratio, the ratio of the 
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number of % at a given pr t to the number of partons(g) at the same pj t : 



7T U 



(j>T t ) 



b n ~ 3 



(8.3) 



I used Eq. 8.3 together with a half and half mixture of quark and gluon jets with b q = 8.2 b g = 11.7 
from Fig. 24b to obtain the parton cross section from the 7T° cross section on Fig. 22b [1 17]. 

8.2.2 Two-particle correlations in A+A collisions are much more complicated than for p-p 

In analogy to Fig. 18 (above), the di-hadron correlations in Au+Au collisions (Fig. 25) show 
a di-jet structure. However, one of the many interesting new features in Au+Au collisions is that 
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Figure 25: a-h) (left) PHENIX [118] azimuthal correlation conditional yield of associated h ± particles 
with p\ for trigger h ± with pf for the various pj p\ combinations shown, i) (right)-(top) PHENIX [119] 
azimuthal correlation function C(A(j>) of h ± with 1 < pj a < 2.5 GeV/c with respect to a trigger h ± with 2.5 < 
Pr t < 4 GeV/c in Au+Au central collisions, where the line with data points indicates C(A<j>) before correction 
for the azimuthally modulated (V2) background, and the other line is the V2 correction which is subtracted 
to give the jet correlation function J(A(j>) (data points), j) (right) -(bottom) PHENIX D parameters [119], 
the angular distance of the apparently displaced peak of the J(A<j>) distribution from the angle A(j) = tc as 
a function of centrality, represented as the number of participants Af part , for the systems and cm. energies 
indicated. 

the away side azimuthal jet-like correlations (Fig. 25c) are much wider than in p-p collisions and 
show a two-lobed structure ("the shoulder" (SR)) at lower pj t with a dip at 180°, reverting to 
the more conventional structure of a peak at 180° ("the head" (HR)) for larger p Tr The wide 
away-side correlation in central Au+Au collisions is further complicated by the large multiparticle 
background which is modulated in azimuth by the V2 collective flow of a comparable width to the 
jet correlation (Fig. 25i). After the V2 correction, the double peak structure ~ ±1 radian from 7T, 
with a dip at K radians, becomes evident. The double-peak structure may indicate a reaction of the 
medium to a passing parton in analogy to a "sonic boom" or the wake of a boat and is under active 
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study both theoretically [120] and experimentally. PHENIX characterizes this effect by the half- 
width D (~ 1 . 1 radian) of the Jet function, /(A0), the angular distance of the displaced peak of the 
distribution from the angle A0 = n. One of the striking features of the wide away side correlation 
is that the width D (Fig. 25j) does not depend on centrality, angle to the reaction plane, pr a and 
y/sNN, which seems problematic to me if the effect is due to a reaction to the medium. 

As discussed above, the xe (also denoted zt) distribution from di-hadron collisions does not 
measure the fragmentation function. Instead, it is sensitive to the ratio of the away parton transverse 
momentum to the trigger parton transverse momentum, xh = pr a / Pr t , which is a measure of the 
differential energy loss of the away parton relative to the trigger parton which is surface biased due 
to the steeply falling p Tt spectrum [121, 122]. The energy loss of the away-parton is indicated by 
the fact that the xe distribution in Au+Au central collisions e.g. for 4 < pj t < 5 GeV/c (Fig. 26a) 
from the data in Fig. 25 [118] is steeper than that from p-p collisions. Eq. 8.2 provides excellent 
fits to these distributions in both p-p and Au+Au collisions (Fig. 26a), where fits to the full away 
side azimuth Head + Shoulder (HS) region are shown. In p-p collisions, the imbalance of the away- 
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Figure 26: a) xe distribution [118] dP/dxE of associated h ± per trigger h ± with 4 < pj t < 5 GeV/c for the 
full away side azimuth, |A0 — n\ < n/2, in p-p (circles) and Au+Au central (0-20%) collisions (squares). 
Fits to Eq. 8.2 are shown with best-fit parameters indicated, b) Same plot for 5 < pr t < 10 GeV/c data [118], 
for the Head region only (HO), \A(j> — 7t\ < 71 / 6. 

parton and the trigger parton indicated by the fitted value of Xh = 0.93 ± 0.03 in Fig. 26a is caused 
by ^-smearing. In A+A collisions, the fitted value x^ = 0.52 ± 0.03 indicates that the away parton 
has lost energy relative to the trigger parton. 

For larger pr t , another interesting effect takes place, punch-through of narrow away-side peaks 
consistent with roughly the same azimuthal width and xe distributions as in p-p collisions [122]. 
Fig. 26b shows the xe distributions and fits at higher 5 < pj t < 10 GeV/c where the region of 
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integration of the away-side azimuthal distribution has been restricted to the head-region only (HO) 
in both p-p and Au+Au collisions. The fitted values of Xh are nearly identical in Fig. 26b (xh = 
0.94 ± 0.01 in p-p and xh = 0.50 ± 0.05) and Fig. 26a but the integrated away-side multiplicity N is 
reduced in the HO region compared to the full-away side (HS). Of more importance to note is that 
the fit to the Au+Au HO data (Fig. 26b) is greatly improved (A^ 2 = 6. 2/1) if a second component 
with the same Xh as the p-p distribution is added (dashed-curve), statistically indicating a parton 
that has apparently punched through the medium without losing energy. The two-component fit 
shows 27% punch-through and a much larger energy loss, 77%. This is more clearly exhibited by 
taking the ratio, Iaa of the xe distributions for Au+Au to p-p (Fig. 27a), and is evident directly in 
the STAR z^(same as xe) distribution [123] (Fig. 27b) by the sharp change in slope at zt (xe) = 0.5 
in Au+Au for 6 < p Tt < 10 GeV/c. 
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Figure 27: a) Iaa(xe) from the ratio of Au+Au/pp distributions in Fig. 26b. Dashed line is 27% punch- 
through from fit. b) STAR [123] away side zt (same as xe) distributions for various pj t in Au+Au, d+Au 
and their ratio in aAaw = 200 GeV collisions. 



The punch-through and normal fragmentation of partons which have lost energy are standard 
features of the QCD energy loss models. For instance the ZOWW [124] model breaks down the 
away-side fragmentation function to: 



D h/C (z c , AE C , M 2 ) = (l-,-<VA>) 



L\ Z 



A/ Zr 



(1) 



where z' c = Pt/(Pt c ~ A£ c ), z' g = (L/A)p T /kE c are the 
rescaled momentum fractions, AE C is the average radiative 
parton energy loss and (L/ A) is the number of scatterings. 
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In words, the three terms in the above equation correspond to normal fragmentation of partons 
which have lost energy; the fragmentation of the lost energy, assumed to be a gluon; normal frag- 
mentation with full original energy of punch-through partons. This 'paradigm' is consistent, as far 
as I can tell, with all measured Iaa{*e) distributions in di-hadron correlations (Fig. 27)-an expo- 
nential decrease at low xe from energy-loss and a flat distribution at large xe due to punch-through 
when no radiation occurs. 

The discrete radiation of relatively large energy gluons is characteristic of the LPM effect [49, 
125] in QCD, where in distinction to QED, it is the radiated gluon that interacts coherently with the 
medium. This suppresses the gluon spectrum for large radiated energy k (rather than at low k as in 
QED, Fig. 4a) but the coherence comes from gluons with a k larger than some critical value ((Qbh) 
which is relatively hard, 1-4 GeV. Thus if LPM radiation occurs there is a large radiation loss, but 
there is a Poisson probability e~ L ^ that no radiation occurs leading to punch-through. However, 
as I understand it, multiple scattering (p T broadening), elastic and normal Bethe-Heitler radiative 
loss should take place even with punch-through, so we still have a long way to go on this issue. 

The energy-loss models all do quite well in describing the suppression, Raa(pt), for single- 
inclusive 7T , although each model has different parameters [126]. However, the suppression of 
direct-e^ from heavy quark decays (recall Fig. 23) disfavors radiative models since, naively, heavy 
quarks should radiate much less than light quarks in the medium. This has attracted much theoret- 
ical attention and is still not explained [127, 128]. 

For two-particle correlations, no clear paradigm has emerged for the two-lobed wide away -jet 
structure. Some suggestions are: 

• Mach or Cerenkov cone due to medium reaction. Since the effect vanishes for larger pj t , with 
no dependence on centrality or pj a (recall Fig. 25j), I find this explanation unlikely. Further 
studies with particle identification might show whether the composition of the shoulder is 
the same as the medium or the jet for a given /7^-kind of like the difference of being hit by 
the wake (water) or the boat (wood). 

• Fluctuations of V3 [129, 130]: (V3) =0 at mid-rapidity but (V3) 7^ 0. This gives the best fits 
to the azimuthal correlation functions C(A0) and the constant value of D « 60° (Fig. 25) and 
is consistent with vanishing at large pj t . Lots of ongoing work at the present time. 

• NLO 3-jet events less suppressed than di-jets due to smaller path-length in medium [131]. 

9. Direct- y-h correlations do measure the fragmentation function 

Direct single- 7 production via the inverse QCD-compton process [132]: 

g + q^Y+q (9.1) 

is an important probe in p-p collisions because it is sensitive the gluon structure function at leading 
order as well as higher orders [133]. It is possibly even more important in A+A collisions because 
the 7 is a direct participant in the hard- scattering (at the constituent level), which emerges from the 
medium without interacting and can be measured precisely. This means that, modulo kj, the pj of 
the quark jet which is equal and opposite to that of the direct-7is known precisely. Furthermore, 
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since the direct-y is not a fragment of a jet, but a direct participant in the hard-scattering, this 
means that the fragmentation function can be measured by the direct- y-hadron xe (zt) distribution, 
as originally suggested by Wang, Huang and Sarcevic in 1996 [134], 'the golden channel' [135]. 

There are nice reviews of both the theory [105] and the experimental discovery of direct- single- 
y production at the CERN-ISR in 1978-80 and other early measurements [136]. Measurements 
of direct-y at RHIC in both p-p [137] and Au+Au [138] collisions have shown the absence of 
suppression for this probe which does not interact with the medium (recall Fig. 23). Calculations 
of direct-y production in QCD give impressive agreement with the data over a wide range of cm. 
energies [133]. Other nice features of direct-y production from the inverse QCD-compton effect 
(Eq. 9.1) is that the y are isolated and the away -jet is a quark-jet, 8 times more likely u than d. 

Beyond Leading Order, single-y may be produced in the fragmentation process [105]; but 
these y are not isolated. However, measurements of y— h azimuthal correlations in p-p collisions at 
both the CERN-ISR [139] and RHIC [140] (Fig. 28) show that direct-y are isolated, with very few, 
if any, accompanying same-side particles, while 7T° have accompanying particles since they are 
fragments of jets from high pr partons. This is strong evidence against a significant fragmentation 
component for direct-y production, easily <10%. 
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Figure 28: Direct single- y-h and 7T°-h azimuthal correlations in p-p collisions: sl)s/snn = 62.4 GeV, pr t > 6 
GeV/c [139]; b) yfs = 200 GeV, p Tt , p Ta as indicated [140]. 

The xe and p out distributions in p-p collisions from 7T°-h and direct- y-h correlations from new 
data at RHIC [140] are shown in Fig. 29. For these measurements of direct- y-h correlations, an 
isolation cut has been made to improve the systematic uncertainty in subtracting the background 
from 7T°-h correlations. The xe distributions show a flatter slope of the 7T°-h distribution compared 
to the direct- y-h distribution, as discussed for Fig. 18d (which is for jet-h vs 7T°-h correlations) [90], 
showing that the y plays the expected role of a jet for correlation measurements, albeit with a much 
higher precision measurement of pj and a much lower rate. Of course, the xe distribution from the 
isolated-direct- y-h measurement is actually the fragmentation function. The p out distributions for 
7T°-h and isolated-direct- y-h correlations (Fig. 29) appear to be nearly identical and both show the 
features of a Gaussian distribution for p out < 3 GeV/c, which is the kj effect, now thought to be due 
to resummation of soft gluons [133], and a power-law tail from NLO hard-gluon emission. This 
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Figure 29: a) xe and b) p oui distributions from 7T°-h and isolated direct- 7-h correlations [140] 

is the first time that the kj effect has been measured with direct- 7; and the yj (k\) is essentially 
identical in 7T° and direct- 7 production. 

9.1 The fragmentation function in p-p collisions 

From the discussion above, the xe (zt) distribution for direct- 7-h (Fig. 29a) should be a mea- 
surement of the u quark fragmentation function. To make this more apparent, we follow the ap- 
proach of Borghini and Wiedemann (Fig. 30a) [141] who proposed using the hump-backed or 
£, = ln(l/z) distribution of jet fragments, which is a signature of QCD coherence [142] for small 




Figure 30: a) Predicted t; distributions in vacuum and in medium for two jet energies [141], together with 
measurements from e+e~ collisions [143, 115] b) t; = lnl/x^ distributions for PHENIX isolated direct-7 
data for all pj t ranges combined, compared to e + e~ collisions at ^Js = 14 and 44 GeV [143]. 

values of particle momentum fraction, z = p/E^ u to explore the medium-modification of jets in 
heavy ion collisions. The use of the t; variable would emphasize the increase in the emission of 
fragments at small z due to the medium induced depletion of the number of fragments at large z. 
The jet energy must be known for this measurement so that it was presumed that full jet reconstruc- 
tion would be required. However, the direct- 7-h correlation is (apart from the low rate) actually 
better for this measurement since both the energy and identity of the jet (8/1 w-quark, maybe 8/2 if 
the q + q 7+ g channel is included) are known to high precision. The PHENIX xe distributions 
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(Fig. 29a) converted to the t, distributions are shown in Fig. 30b in quite excellent agreement with 
the dominant w-quark fragmentation functions measured in e + e~ collisions at \/~s/2 = 7 and 22 
GeV [143], which cover a comparable range in jet energy. Also, Fig. 30b is plotted on a log scale 
for dN/dt; which allows the full range of the fragmentation function, notably for z > 0.2, £ < 1.6 
to be visible. 



9.2 "The golden channel", direct- y-h in Au+Au collisions 

Before discussing direct- y-h correlations in Au+Au collisions, it is worth reviewing the re- 
sults for 7T°-h correlations, presented in Fig. 31 [144] as Iaa, the ratio of the pr distributions of 
the partner-h in Au+Au central collisions to p-p collisions, for several values of trigger 7F° pj t at 
y/sNN = 200 GeV. All the distributions show the characteristic shape of an exponentially falling Iaa 
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Figure 31: a) PHENIX [144] Away-side 7 A a for a "head-only" (HO) \A<j> - n\ < n/6 selection (solid 
squares) and the entire away-side (H+S), |A0 — 7t\ <7t/2 (solid circles) as a function of partner-h pr for var- 
ious pr t . Calculations from two different predictions are shown for the head region in applicable partner-h 
Pt ranges. Raa(pt) of 7T° for pj t > 5 GeV/c [145] are included for comparison. 

for lower values of pj a due to energy loss of the away -jet and a flat distribution for larger values of 
pr a due to partons that punch-through the medium with no energy loss and then fragment normally 
as in p-p collisions. Two other important points are evident from Fig. 31: a) the H+S and HO data 
become the same for pj t > 1 GeV/c because the away-jets in Au+Au are no longer anomalously 
wide and are consistent with the shape of the p-p azimuthal correlation in both the falling and flat 
part of the distribution; b) In the punch-thru region, Iaa > Raa at the same value of p?. 

The expectations for the difference in behavior of 7T°-h and direct- y-h correlations for the cases 
of either a totally absorbing medium (surface emission only) and a translucent lossy medium (par- 
tons can emerge from the medium having lost energy) are briefly summarized: 

Totally absorbing medium: 

direct-y-h: y emerge from all throughout the medium. Away partons are either totally absorbed 
or emitted from the surface with normal fragmentation, same surface bias as for inclusive 
hadrons (but on the away surface). Iaa(pt^) =constant=fraction not absorbed^ R\ A (pr t ). 
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7T-h: For inclusive hadrons, strong surface bias: only partons emitted from the surface are seen. 
For di-hadrons, only tangential emission is seen. Iaa =constant=the fraction of surface emis- 
sion that is tangent to the surface. Therefore Iaa < ^aa(^)- 

Translucent lossy medium: 

direct- y-h: Away partons lose energy. Some punch through or are tangent. Iaa is exponentially 
decreasing at small xe due to partons with energy loss, constant at larger xe due to punch- 
thru. 

7T°-h: Roughly the same except away partons may lose more energy due to surface bias of the 
trigger, and fewer may punch through. 

The ;r -h data in Fig. 3 1 show Iaa > Raa so disfavor the totally absorbing medium. The present 
direct- y-h data are less conclusive. 

Fig. 32a shows the PHENIX preliminary xe (zt) distributions for direct- y-h correlations in 
p-p and Au+Au from QM2009 [146]. The p-p data show nice xe scaling for all pr t (as would 
be expected for a fragmentation function) and are fit to an exponential, e~ bzT , with slope b = 
6.89 ±0.64. The Au+Au data also exhibit reasonable xe scaling for all pj t with an exponential 
slope b = 9.49 ± 1.37. The xe scaling and steeper slope in Au+Au suggest a constant fractional 
energy loss of the away-side parton; but the statistical significance is marginal. The thin red line 
shows that possibly b —> 6.89 for zt > 0.4 which would imply /^A^constant for xe > 0.4 as for 
7r°-h. The STAR [147] measurement of 7T -h and direct- y-h Iaa as a function of zt for 8 < pj t < 16 
GeV/c are shown in Fig. 32b and show equal and constant Iaa — 0.3 for 0.3 < zt < 0.9 for both 
the 7T°-h and direct-y correlations. There is no evidence in either case for an exponential rise of Iaa 
towards zt = 0, possibly because the zt > 0.3 range is too high. Also the STAR value of Iaa — 0.3 
for 7T°-h is smaller than the PHENIX value, Iaa = 0.50 ± 0.08, in the same pr t range. 



a) 




Figure 32: a) PHENIX [146] preliminary zt fe) distribution (|A0 — 7t\ < n/6) for direct- y-h, in Au+Au 
(red), and multiplied by a factor of 10 for p-p (blue), at ^snn = 200 GeV. b) STAR [147] away-side Iaa 
vs. zt for 71° -h (triangles) and direct- y-h (circles) with 8 < pr t < 16 GeV/c at at ^s^n = 200 GeV. For a 
discussion of the curves see Refs. [147, 135]. 
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In summary, the results of measurements of the "golden" (but difficult) channel, direct- y-h, are 
at too early a stage to tell whether this channel will live up to its promise to understand the QGP, 
and QCD in a medium; i.e. whether it will be a lasting paradigm, or go the way of the J/yf. 
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